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After a brief summary of the QCD phase structure with light quarks, we discuss two re- 
cent developments on in-medium hadrons. First topic is the cr-meson which is a fluctuation 
of the chiral order parameter qq. Although a is at best a broad resonance in the vacuum, 
it may suffer a substantial red-shift and show a characteristic spectral enhancement at the 
2m n threshold at finite temperature and baryon density. Possible experimental signatures 
of this phenomenon are also discussed. Another topic is the first principle lattice QCD cal- 
culation of the hadronic spectral functions using the maximum entropy method (MEM). 
The basic idea and a successful example of MEM are presented. Possible applications of 
MEM to study the in-medium hadrons in lattice QCD simulations are discussed. 

1. Introduction 

One of the most intriguing phenomena in quantum chromo dynamics (QCD) is the 
dynamical breaking of chiral symmetry. This explains the existence of the pion and 
dictates most of the low energy phenomena in hadron physics. The dynamical breaking 
of chiral symmetry is associated with the condensation of quark - anti-quark pairs in 
the QCD vacuum, (qq), which is analogous to the condensation of Cooper pairs in the 
theory of superconductivity As the temperature (T) and/or the baryon density (rig) 
increase, the QCD vacuum undergoes a phase transition to the chirally symmetric phase 
where (qq) vanishes. Studying what are the possible phase structure and how the phase 
transition takes place are one of the main aims of the modern hadron physics ||. 

When we talk about the phase structure, we need to fix appropriate variables to define 
the phases. Such variables in QCD are current masses of light quarks {m u ^ s ), T, and rig. 
Since the following relation holds, neglecting m u ^ in the first approximation and studying 
the phase structure in the m s -T-riB space would give us a good insight into the real world: 

m u ,d < ™ s ~ T c ~ n l J 3 ~ A QCD - (1) 



Here T c (n c ) are the critical temperature (density). In Fig.l, possible phase structures in 
the m s -T plane (with ne=0) and m s -nB plane (with T=0) are shown. m s =0 (m s =oo) 
corresponds to the limit of SU(3)®SU(3) (SU(2)®SU(2)) chiral symmetry. Recent lattice 
simulations show T c ~ 175 (155) MeV for 3 (2) flavors [[J. Furthermore, the first order 
transition for small m s turns into the second order one for large m s at the tricritical point 
[|J. The precise value of m s at the tricritical point is not known yet. 
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Figure 1. Possible phases in the m s -T plane at ub=0 (left panel) and the m s -nB plane 
at T=0. m u ^=0 is assumed. Relevant degrees of freedom in each phase are also shown. 
uq = 0.17fm -3 is the normal nuclear matter density. 



The phase transition at finite baryon density is not well understood mostly because 
the information from lattice QCD is poor. Naive analysis based on the bag model shows 
the first order transition ||, which implies the existence of a mixed phase in the m s -nB 
plane as shown in Fig.l (right panel). However, the possibility of a smooth hadron-quark 
transition is not excluded. The bag model analysis also suggests the strange quark matter 
as the true ground state of matter for sufficiently small m s || . The system at finite baryon 
density, which is intrinsically quantum, is obviously richer in physics than that at finite 
T and has various phases such as the color superconductivity, meson condensations, and 
baryon superfluidity (see the reviews, 0). 

Now, what are the observables associated with the chiral transition? Model independent 
statement is that hadronic correlations in a same chiral multiplet (such as S (scalar)— P a 
(pseudo-scalar) and (vector)— A® (axial- vector)) should be degenerate when (qq) — > 0, 
namely 

(S(x)S(y)) -> (P a (x)P a (y)), (A%x)A>(y)) -> (V«(x)V u \y)). (2) 
Thermal susceptibilities for hadronic operators O defined in the Euclidian space, 

rl/T r 

X H =J Q drj d 3 x(O\r,x)O(0,0)), (3) 

may be used to check such chiral degeneracy. Shown in Fig. 2 is the full QCD simulation 
of \Jl/x H on ^ ne lattice with 2 light flavors 0. One can see the degeneracy between 
the cr-channel (1=0, J p =0 + ) and the 7r-channel (1=1, J p =0~) above the critical point. 
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Figure 2. Thermal susceptibilities in three different channels (tt, a, and 5) for two-flavor 
QCD on the 8 3 x 4 lattice with m u ^ a = 0.02 ||. The vertical (horizontal) axis denotes 

\Ji/Xh (the lattice coupling (3 = Q/g 2 ). 



There remains a splitting between the 5-channel {1=1, J p =0 + ) and the 7r-channel above 
the critical point, which reflects the breaking of Ua{1) symmetry. 

To make a direct connection with the experimental observables, however, one needs to 
go further and calculate the Minkowski (real-time) correlation and the hadronic spectral 
functions at finite T and @. The aim of this talk is to show some recent attempts 
toward this goal. In Sec. 2 and Sec. 3, the spectral function in the cr-channel will be 
discussed near 2n threshold where almost model-independent argument is possible. 

Experimentally, the resonances in the scalar and axial-vector channels have large width 
and are difficult to measure with possible exception discussed in Sec. 2, 3. For this reason, 
in-medium neutral vector-mesons which couple directly to the virtual photon have been 



studied both theoretically and experimentally as a probe of the chiral transition |L0 
Unfortunately, theoretical calculations of the spectral function in the vector channel on 
the basis of the effective Lagrangians are rather model dependent, and the first principle 
lattice QCD study has been called for. Recently, a new promising approach was proposed: 
it utilizes the maximum entropy method for extracting the spectral functions from lattice 
QCD data. We will discuss this new development in Sec. 4. 

2. Spectral enhancement in hot plasma 

The fluctuation of the order parameter becomes large as the system approaches to 
the critical point of a second order or weak first order phase transition. In QCD, the 
fluctuation of the phase (amplitude) of the chiral order parameter (qq) corresponds to 
7r (cr). The vital role of such fluctuations in the dynamical phenomena near the critical 
point at finite T was originally studied in |p"T| , |T2"| . 



It was suggested that the chiral restoration gives rise to a softening (the red-shift) of 
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Figure 3. Spectral functions in the 7r channel (A) and in the a channel (B) for T=50, 
120, and 145 MeV. |15] p a (uS) in (B) shows only a broad bump at low T (a), while the 
spectral concentration is developed as T increases, (a)— >(b)— >(c). 



the a, which in turn leads to a small a- width due to the suppression of the phase-space 
of the decay o 



2tt 11 



Therefore, a may appear as a narrow resonance at finite T 
although it is at best a very broad resonance in the free space with a width comparable 
to its mass [jn|. (For the phenomenological applications of this idea in relation to the 
relativistic heavy ion collisions, see |]T4"|.) 

Further theoretical analysis by taking account the coupling a <-> 2n shows that (i) the 
spectral function is the most relevant quantity for studying the true nature of a, and 
(ii) the spectral function of o has a characteristic enhancement at the 2n threshold near 
T c |15||. In Fig. 3, shown are the spectral functions p 7r ( cr ) in the n (cr)-channel at finite T 
calculated in the 0(4) linear a model: Two characteristic features are the broadening of 
the pion peak (Fig. 3 (A)) and the spectral concentration at the 2n threshold (u ~ 2m n ) 
in the cx-channel (Fig.3(B)). The latter may be measured by the 27 spectrum from the 
hot plasma created in the relativistic heavy ion collisions Jl5] . 



3. Spectral enhancement in nuclear matter 

(qq) at finite baryon density obeys an exact theorem in QCD [jn 



(gg) 
(qq)o 



n B 



T, nN + in— 
dm 



d ( EOhb) 



A 



(4) 



where H n N = 45 ± 10 MeV is the pion-nucleon sigma term and E(ub)/A is the nuclear 
binding energy per particle with m = (m u + m<f)/2. (qq)o denotes the chiral condensate 
in the vacuum. The density-expansion of the right hand side of (|j) gives a reduction of 
almost 35 % of (qq) already at the nuclear matter density n = 0.17fm~ 3 . 
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Figure 4. In-medium tttt cross section in the I—J—Q channel as a function of the cm. 
energy yfs for several different values of the condensate ratio $ = (a)/(a) in nuclear 
matter 



The near-threshold enhancement discussed in Sec. 3 has been also studied at finite 
baryon density |17j]. Because of the decrease of the chiral condensate discussed above, 
the spectral function in the cr-channel, p a , has a following generic behavior at densities 
not far from no; 



p a (u ~ 2m n ) oc 



6(u - 2m n ) 



„,2 



(5) 



which shows a spectral concentration at the 2n threshold. 

Recently CHAOS collaboration [19] reported the data on the tt + tt ± invariant mass 
distribution M^ +7r± in the reaction A(tt + , tt + tt ± )A' with the mass number A ranging from 
2 to 208: They observed that the yield for M£+ n _ near the 2tt threshold is close to zero 
for A = 2, but increases dramatically with increasing A. They identified that the tt + tt~ 
pairs in this range of M^ +7r _ is in the I=J=0 state. This experiment was originally 
motivated by a possibility of strong tttt correlations in nuclear matter [fJU] . However, the 



state-of-the-art calculations using the nonlinear chiral Lagrangian together with 7riVA 
many-body dynamics do not reproduce the cross sections in the 1=0 and 1=2 channels 
simultaneously [^T|; the final state interactions of the emitted two pions in nuclei give rise 
to a slight enhancement of the cross section in the 1=0 channel, but is not sufficient to 
reproduce the experimental data. This indicates that an additional mechanism such as 
the partial restoration of chiral symmetry may be relevant for explaining the data |TT|,|2"2" 



To make a close connection between the idea of the spectral enhancement and the 
CHAOS data, the in-medium tt-tt cross section has been calculated in the linear and non- 



linear a models ||23|| . It was shown that, in both cases, substantial enhancement of the 
tt-tt correlation in the I=J=0 channel near the threshold can be seen due to the partial 
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restoration of chiral symmetry. Furthermore, an effective Air-N-N vertex responsible for 
the enhancement in the non-linear chiral Lagrangian is identified. In Fig. 4, the in-medium 
7T-7T cross section in the cx-channel calculated in the 0(4) linear a model is shown. 

Theoretically, it is of great importance to make an extensive analysis of the n-n in- 
teraction in heavy nuclei with the effect of partial chiral restoration to understand the 
CHAOS data. Experimentally, it is definitely necessary to confirm the CHAOS result not 
only in the same reaction but also in other reactions. Measuring a — > 2tt° — ► 47 is one of 
the clean experiments, since it is free from the p meson background inherent in the 7r + 7r~ 
measurement. Measuring a — > 27 is interesting because of the small final state interac- 
tions. Dilepton detection through the scalar-vector mixing in matter, a — > 7* — > e + e~ 
and the formation of the a mesic nuclei through the nuclear reactions such as (d, t) are 



other possible experiments [ITS jl? 



4. Spectral functions from lattice QCD 

The lattice QCD simulations have remarkable progress in recent years for calculating the 
properties of hadrons as well as the properties of QCD phase transition In particular, 
masses of light mesons and baryons in the quenched QCD simulation agree within 5-10 % 
with the experimental spectra |]25| . However, the lattice QCD had difficulties in accessing 
the dynamical quantities in the Minkowski space, because measurements on the lattice 
can only be carried out for discrete points in imaginary time. The analytic continuation 
from the imaginary time to the real time using the noisy lattice data is highly non-trivial 
and is even classified as an ill-posed problem. 

Recently a new approach to extract spectral functions of hadrons from lattice QCD 
data by using the maximum entropy method (MEM) was started ||26|| . MEM has been 



successfully applied for similar problems in image reconstruction in crystallography and 
astrophysics, and in quantum Monte Carlo simulations in condensed matter physics [p? 



There are three important aspects of MEM: (i) it does not require a priori assumptions 
or parametrizations of SPFs, (ii) for given data, a unique solution is obtained if it exists, 
and (iii) the statistical significance of the solution can be quantitatively analyzed. 

4.1. Basic idea of MEM 

The Euclidean correlation function D(t) of an operator 0(t, x) and its spectral decom- 
position at zero three-momentum read 

T,x)O(0,0))d 3 x = J K(T,u)A(Lu)du, (6) 

where u is a real frequency, and A(u) is the spectral function (or sometimes called the 
image), which is positive semi-definite. K(t,u) is a known integral kernel (it reduces to 
Laplace kernel e~ TU} at zero T.) 

Monte Carlo simulation provides D{ji) on the discrete set of temporal points < Tj/a < 
A^. From this data with statistical noise, we need to reconstruct the spectral function 
A(u) with continuous variable u. This is a typical ill-posed problem, where the number 
of data is much smaller than the number of degrees of freedom to be reconstructed. 
This makes the standard likelihood analysis and its variants inapplicable unless strong 
assumptions on the spectral shape are made [p8fl . MEM is a method to circumvent this 



7 



difficulty through Bayesian statistical inference of the most probable image together with 
its reliability. 

Let's start with the Bayes' theorem: P[X\Y] = P[Y\X]P[X]/P[Y], where P[X\Y] is 
the conditional probability of X given Y. The most probable image A(u) for given lattice 
data D is obtained by maximizing the conditional probability 



P[A\D) oc P[D\A]P[A\. 



(7) 



The reliability of the obtained result can be checked by the second variation of PL4|.D] 
with respect to A(uj). 

P\D\A] in (0) is a standard x 2 ? namely P[D\AH] = Z^ 1 exp(-L) with a likelihood 
function L = (1/2) J2ij(D( T i) — D A {ji)) C^ 1 (D{jj) — D a (tj)) and a normalization fac- 
tor Z L . D{ji) is the lattice data averaged over gauge configurations and P A (Tj) is the 
correlation function obtained by given A. C is a covariance matrix of the data. 

P[P|A] (the prior probability) in (|7p, can be written with parameters a and m as 
P[A) = Zg 1 exp(aS') by the combinatorial argument. Here S is the generalized informa- 
tion entropy, 



S 



A{u) 



m[uj) 



A{uu) log 



■A( U ) 



m w 



(8) 



with Zs being a normalization factor, a is a real and positive parameter and m(uj) is a 
real and positive function called the default model. 

The output image A out is given by a weighted average over A and a: 



A 



out 



A a {ui) P[a\Dm] da, 



(9) 



where A a {uj) is obtained by minimizing the "free-energy" F = L — aS. a dictates the 
relative weight of the entropy S (which tends to fit A to the default model m) and the 
likelihood function L (which tends to fit A to the lattice data). Note that a appears only 
in the intermediate step and is integrated out in the final result. 

One can prove that the solution of 5F = is unique if it exits. The error analysis of the 
reconstructed image can be studied by evaluating the second derivative of the free energy 
5 2 F/5A(lj)5A(uj'). The default model m is chosen such that the error becomes minimum. 

4.2. MEM with lattice data 

In Fig. 5, shown is an example of the spectral function in the pion and the rho meson 
channels at T = extracted from the quenched lattice QCD data [26[ . The lattice size is 
20 3 x 24 with (3 = 6.0, which corresponds to a = 0.0847 fm (a" 1 = 2.33 GeV), and the 
spatial size of the lattice L s a = 1.69 fm. Hopping parameters are chosen to be k = 0.153, 
0.1545, and 0.1557 with N con f = 161 for each k. For the quark propagator, the Dirichlet 
(periodic) boundary condition is employed for the temporal (spatial) direction. We use 
data at 1 < Tj/a < 12. 

The obtained images have the low-energy peaks corresponding to 7r and p, and the broad 
structure in the high-energy region. The mass of the p-meson in the chiral limit extracted 
from the peak in Fig. 5 is consistent with that determined by the asymptotic behavior of 
D(t). Although the maximum value of the fitting range T max /a = 12 marginally covers 
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Figure 5. Spectral functions p out {uj) = A out (u>)/u> 2 obtained by MEM using the quenched 
lattice QCD data. The lattice size is 20 3 x24 with a = 0.0847 fm. 12 data points in the 
temporal direction (1 < Tj/a < 12) are used for the MEM analysis, (a) is for the pion 
channel and (b) is for the rho-meson channel. The figures are taken from |21| . 



the asymptotic limit in r, we can extract reasonable masses for ir and p. The width of re 
and p in Fig. 5 is an artifact due to the statistical errors of the lattice data. In fact, in the 
quenched approximation, there is no room for the p-meson to decay into two pions. 

As for the second peaks, the error analysis shows that their spectral "shape" does not 
have much statistical significance, although the existence of the non-vanishing spectral 
strength is significant. Under this reservation, we fit the position of the second peaks and 
made linear extrapolation to the chiral limit and obtain results consistent with experi- 
mental data on ir' and p'. 

The error in Fig. 5 indicates the uncertainty of the spectrum averaged over the interval 
in the frequency space. Better data with smaller a and larger lattice volume will be helpful 
for obtaining spectral functions with smaller errors. 



4.3. The future of MEM 

MEM introduces a new way of extracting physical information from the lattice QCD 
data. It is now time to study how hadrons are modified in the medium using this new 
approach. The long-standing problem of in-medium spectral functions of vector mesons (p, 
u>, 0, J/ip, T, • • •, etc) and scalar/pseudo-scalar mesons (ex, vr, • • •, etc) can be studied using 
MEM combined with finite T lattice simulations. The in-medium behavior of the light 
vector mesons [21]] and scalar mesons [[□]] is intimately related to the chiral restoration in 
hot/dense matter, while that of the heavy vector mesons is related to deconfinement [|30[| . 
Simulations with an anisotropic lattice is necessary for this purpose to have enough data 
points in the temporal direction at finite T, which is now under way [|3T|. (See also an 
attempt on the basis of NRQCD simulation p2| .) It is interesting to study the spectral 
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functions with finite three-momentum k for treating moving-hadrons in the medium. 
Non-perturbative collective modes above the critical temperature of the QCD phase 



transition speculated in |l ^|3^j34[ may be studied efficiently by using MEM, since the 
method does not require any specific ansatze for the spectral shape. Correlations in the 
diquark channels in the vacuum and in medium are interesting to be explored in relation 
to the q-q correlation in baryon spectroscopy and to color superconductivity at high Ub 



5. Summary 

Although the light a-meson does not show up clearly in the free space because of its 
large width due to the strong coupling with two pions, it may appear as a soft and narrow 
collective mode in the hadronic medium when the chiral symmetry is (partially) restored. 
A characteristic signal is the enhancement of the spectral function and the it-it cross 
section near the 2tt threshold. They could be observed in the hadronic reactions with 
heavy nuclear targets as well as in the heavy ion collisions through the detection of pion 
pairs and photon pairs. Finding such signal provides us with a better understanding of 
the non-perturbative structure of the QCD vacuum and its quantum fluctuations. 

The maximum entropy method (MEM) opens the door to extract the spectral func- 
tions from the first principle lattice QCD data without making a priori assumptions on 
the spectral shape. The method has been tested at T=0 with reasonable success. The 
uniqueness of the solution and the quantitative error analysis make MEM superior to any 
other approaches adopted previously. The method has a promising future for analyzing 
in-medium properties of hadrons and its relation to chiral structure of matter. 
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